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THE DETERMINISTIC LIMIT OF THE MORAN MODEL: 
A UNIFORM CENTRAL LIMIT THEOREM 

FERNANDO CORDERO^ 


Abstract. We consider a Moran model with two allelic types, mutation and 
selection. In this work, we study the behaviour of the proportion of fit indi¬ 
viduals when the size of the population tends to infinity, without any rescaling 
of parameters or time. We first prove that the latter converges, uniformly in 
compacts in probability, to the solution of an ordinary differential equation, 
which is explicitly solved. Next, we study the stability properties of its equi¬ 
librium points. Moreover, we show that the fluctuations of the proportion of 
fit individuals, after a proper normalization, satisfy a uniform central limit 
theorem in [0,oo). As a consequence, we deduce the convergence of the corre¬ 
sponding stationary distributions. 
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Mathematics Subject Classification (2010): Primary 92D25, 60F05; Sec¬ 
ondary 60J27, 60J28 


1. Introduction 

One of the main goals in population genetics is to describe the evolution of a 
population subject to the evolutionary forces of mutation and selection. Significant 
work has been done in this direction (for review see [H El El 13 ) ■ Among a large 
variety of stochastic and deterministic models, finite population models are consid¬ 
ered to be more realistic. However, it is not always easy to deal with them. This 
situation leads to consider large population approximations. 

In this paper we consider a Moran model with population size two allelic 
types, mutation, and selection. Rather than looking to its classical diffusion limit 
approximation, we aim here to study the behaviour of the proportion of fit indi¬ 
viduals when the size of the population tends to infinity, without any rescaling of 
parameters or time. The theory of density dependent families of Markov chains 
provide us with the necessary tools to achieve this goal (see El uni HH El). 

In a first step, as an application of El Theorem 3.1], we show a dynamical 
law of large numbers, which tell us that the proportion of fit individuals converges, 
uniformly in compacts in probability, to the solution of a ordinary differential equa¬ 
tion. The latter is explicitly solved and an analysis of the stability of its equilibrium 
points is given. We also establish the connection of the deterministic limit with the 
well-known deterministic 2-allele parallel mutation-selection model (see [3]). Next, 
we turn out to the study of the fluctuations of the proportion of fit individuals 
around its deterministic limit. More precisely, based on [TUI Theorem 2.7] and [TTl 
Theorem 8.5], we prove that, after a proper normalization, these fluctuations be¬ 
have asymptotically as a Gaussian process. This result takes the form of a uniform 
central limit theorem in [0, oo). In [111 Theorem 8.5], the corresponding Gaussian 
process is characterised as the solution of a stochastic differential equation, which 
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in our setting can be explicitly solved. Finally, as a consequence of the previous 
results, we prove that the underlying stationary distributions converge to a Dirac 
mass at the unique stable equilibrium point of the aforementioned ordinary differ¬ 
ential equation. 

The paper is organized as follows. In Section ^ we shortly describe the 2-type 
Moran model with selection and mutation and we recall its basic properties. In 
Section [21 we state and prove our main results, which deal with the asymptotic be¬ 
haviour of the proportion of ht individuals: (1) a law of large numbers - Proposition 
(Ih.lD . (2) a uniform central limit theorem - Theorem (liI.4H . and (3) the convergence 
of the stationary distributions - Corollary (13.611 . 


2. Some basic facts on the 2-type Moran model 


The 2-type Moran model of size N with selection and mutation describes the 
evolution in (continuous) time of a population of size N in which each individual is 
characterised by a type i G {0,1}. The underlying dynamic is given as follows. If 
an individual reproduces, its single offspring inherits the parent’s type and replaces 
a uniformly chosen individual, possibly its own parent. The replaced individual 
dies, keeping the size of the population constant. 

Individuals of type 1 reproduce at rate 1, whereas individuals of type 0 reproduce 
at rate 1 + sn, sn > 0. Mutation occurs independently of reproduction. An 
individual of type i mutates to type j at rate un Vj, un > 0, vj € (0,1), pq + i^i = 1- 
We denote by the number of individuals of type 0 at time t. The process 
:= (Af) t>o is a continuous-time Markov chain with infinitesimal generator 

Ax-fik) := (/(fc + 1) - /(fc)) + {fik - 1) - /(fc)), (2.1) 


where 


k (1 + Sn) + {N — k) UN Po if ^ = 1, 




— < + kuNVi 


0 


if^ = -l, 
if |€| > 1, 


and := —q^k+i ~ ^kk-i- other words, X^ is a birth-death process with 
birth rates '■= q^k+i death rates ^k ■= ^kk-i- particular, when un > 0, 
X^ admits a unique stationary distribution, which is given by 


TTxHk) 

tl 


ke{o,...,N}. 


where Cn is a normalising constant (see i)- When Un = 0, by contrast, A^ is an 
absorbing Markov chain with 0 and N as absorbing states. 

It is well known that, when the parameters of selection and mutation satisfy 


lim NuN = dG{0,oo) and lim Asjv = cr G (0, oo), (2.2) 

N —¥00 N —¥00 

the rescaled process {X^^/N)t>o converges in distribution to the Wright-Fisher 
diffusion (see, e.g., [3 p. 71, Lemma 5.11]). By the latter we mean the continuous¬ 
time Markov process Y := (Yt)t>o with infinitesimal generator given by 

Ay fix) := a:(l — x)^^j^ix) + (tT(l — x)x + Oi'oil — x) — Ovix) ^{x), x G [0,1]. 
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3. The deterministic limit 

3.1. A law of large numbers. In contrast to the diffusion limit framework, where 
it is assumed that the parameters of the model satisfy (EH), we consider here 
constant parameters of mutation and selection, i.e. un = u > 0 and sat = s > 0. 
In addition, we do not rescale the time. In this setting, a deterministic limit emerges 
when the size of the population converges to infinity. To see this, we first observe 
that the infinitesimal parameters of satisfy 

<lk,k+e = , ^^0, 

where g : R x Z \ {0} —>■ K is defined by 

( (l + s)p(l-p) + Mi/o(l-_p) iff = l, 
q{p,£):=< p{l-p)+ui/ip if£=-l, 

[ 0 if \i\ > 1. 

We conclude that the sequence of Markov chains is density dependent 

in the sense of 0 Sect. 3] (see also [TO])- Let us denote by := {Z^)t>o the 
continuous-time Markov chain given by 

:= > t > 0. 

t N ~ 

In addition, we set F := q{-, I) — q{-, —I), i.e. 

F{x) := sa;(I — x) + Mr'o(I — x) — uvix = —sx^ + (s — u)x + mz/q, a; G M. 

The next result provides the asymptotic behaviour of the process Z^ in the form 
of a dynamical law of large numbers. 

Proposition 3.1 (Law of large numbers). For each zq G [0,1], we denote by z{zo, •) 
the solution of 

dz 

z(0) = zo G [0,1] and —{t)=F{z{t)), tG[0,T]. (3.1) 

Assume that limAr_>.oo Zq = zq € [0,1]. Then, for all e > 0, we have 
lim P (sup\Z^ — z{zo,t)\ > e] = 0, 

M-s-oo \t<T j 

i.e. Z^ converges to the solution of (EH) uniformly in compacts in probability 
(convergence ucp). 

Proof. Note that the function F is Lipschitz in [0,1] and that F{Q) = mi/q > 0 and 
F(l) = —uvi < 0. We conclude that, for all zq G [0,1], Eq. (13.11) has a unique 
solution z{zo, •) defined in [0, oo), which in addition does not leave the interval [0,1]. 
Besides that, q{x,t) = 0 when |I'| > 1, and hence 

sup \i\q{x,i) < oo and lim sup ^ |^|g(a;, ^) = 0. (3.2) 

The desired result follows therefore as an application of Theorem 3.1] (see also 
[SI Theorem 11.2.1]). □ 

Remark 3 . 2 . An analogous result has been shown in [5] for a discrete-time version 
of the Moran model. 

Note that when s = u = 0, F is identically equal to 0. Thus, the solution 
of (13.11) is constant. When s = 0 and u > 0, Eq. (13.11) is linear and has a unique 
equilibrium point given by vq. In the remaining case, Eq. (13.11) is a Riccati equation 
with constant coefficients and can hence be solved by means of quadratures (see 
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m)- Moreover, en) has two equilibrium points, Xg < 0 < Xij" < 1, given by the 
zeros of F, i.e. 


Xn := 


s — u —'/a , s — M + 

and xiT := 


2s 


2s 


(3.3) 


where A := (s — FAsuvq. For simplicity, when s = 0 and u > 0, we set x(j" := vq 
and A := u. 

The following lemma gives the explicit expression of the solution of (EH) and 
the stability properties of its equilibrium points. 


Lemma 3.3. For each zq G [0,1], the solution of (13.11) is given by 
' zo if s = u = 0, 

Vo + [zo — i/q) if s = 0 and m > 0, 

x+{zo-x~)-x~(zQ-x+)e '’(“'o “=0 )* 


z{zo,t) := < 


= (x+-x^ )t 


i/s > 0, 


{zo-x^ )-{zo-xf)e 

for all t > 0. Moreover, when s + it > 0, the equilibrium point Xq is asymptotically 
stable. More precisely, 

lim z(zo,t) = Xq . (3.4) 


In addition, when s > 0, the equilibrium point Xg is unstable. 


Proof. The first part of the statement is obtained through a straightforward veri¬ 
fication. When s + u > 0, eh is obtained from a direct calculation of the corre¬ 
sponding limit, or simply noting that 

dF 

— (x+) = -A<0. (3.5) 

Similarly, if s > 0, ^{xf) = ^fK > 0. The last assertion follows. □ 


3.2. The underlying deterministic mutation-selection model. In this sec¬ 
tion, we shortly describe the well-known deterministic 2-allele parallel mutation- 
selection model and its relation with Eq. EH- This model goes back to H p. 
265] and is usually obtained via direct deterministic modelling, rather than via a 
limiting procedure starting from a stochastic model. 

The deterministic 2-allele parallel mutation-selection model describes the evolu¬ 
tion of a population in which each individual is characterised by a type i G {0,1}. 
The absolute frequencies of individuals of type 0 and 1 satisfy the following system 
of differential equations 

where yl := F + * “ “"' , ). (3.6) 

dt^ ’ ^ " V Wo l-uvoj ^ ' 


It is straightforward to see that, under the initial condition ?/(0) := {zo,l — zq) 
with Zq G [0,1], the proportion of individuals of type 0 evolves following Eq. (13.ip . 
Therefore, we have 

yo{t) 


z{zo,t) = 


(3.7) 


yo{t)+yi(t)' 

Similarly, Eq. EH can be turned into the linear system of differential equations 
EH using the following transformation (see El, m or m) 

y{t) ■= (yoit), yiit)) = exp I t -l- s J z{u)du I X (z(t), 1 - z{t)). 
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We point out that this correspondence between the solutions of 63) and (63 
gives an alternative way to obtain Lemma 63 We provide the details in the case 
s > 0, the other cases are obtained in a similar way. We first express y as 

y{t)=y{Q)e^\ 

Additionally, the eigenvalues of the matrix A are A+ := l + sa;[j" and A_ := l + sa;[^, 
with corresponding right eigenvectors given, up to multiplicative constants, by 


f Wo + sa;Q 

?;+ := 

V UI/Q 


and + 

uvo 


Thus, the matrix can be explicitly computed through diagonalization of the 
matrix A. In particular, with the initial condition j/(0) = ( 20 , 1 — zq), we obtain 

o(l + SX^ 


X+)t 


yo{t) = — -^ [xoizo-X q) - Xq{zo- x^)e 

XA — Xr. \ 


r)‘ 


Xo - X, 
,{l+sx, 


)■ 


yiit) = ^ ((1 - - a;o ) - (1 - Xo ){zo - x^)e ^0 . 

Plugging this in (13.71) yields (13.311 . This alternative way of solving equation (13.11) 
has the advantage that it can be easily generalized to the multi-type case. 

3.3. A uniform central limit theorem in [0, 00 ). We assume in the sequel that 
M > 0. In this setting, it would be natural to study the behaviour of the stationary 
distribution of , denoted by tt^n, when N tends to infinity. The convergence 
ucp of to 2 , derived in Section 63 is not sufficient to deduce such a result. For 
this reason, we aim to provide a uniform convergence result on [0, 00 ) providing the 
convergence of the limiting distributions. The central limit theorem in m (see also 
[mug) for density dependent families of Markov chains give us a way to achieve 
our task. To see this, we first define 

g{x) := q{x, 1) -|- q{x,—l) = —(2 -|- s)x'^ -|- (2 -|- s — u{vo — vi)) x + Wq, x G R. 
Next, we consider, for zq G [0,1], the Gaussian diffusion V^° := {Vf°)t>o given by 

Pi^izoA)) f dBy ifzoT^x;} 


:= 


\j9{xt)e 


F(z(zo,ti)) 
t 

-VAt r 




/e' 

0 


’ dBy \izo = XQ, 


where (Bt)t>o is a standard Brownian motion and z(zo,-) is the solution of (I3.1|l . 
One can easily check that, for all x G [0,1], g{x) > u{vo twi) >0. Hence, V^° is 
well defined. 

Now, we introduce the following characteristic functions: 

'^i9VN{Z^-z(zo,t))' 


-E 


and ip(t, 9) := E 


jsvr 


Theorem 3.4 (Central limit theorem). Assume that limAr_,.oo 
Then \/N{Z^ — z(zo, •)) — > V^°. Moreover, we have 


t>0,9eR. 
y/N{Z^ - zo) = 0. 


N —¥oo 


lim supl'ipNitjO) — tjj{t,9)\ = 0. 

AT-J-oo t>Q 

Proof. We focus only on the case zq Xq . The case zo = Xq is proven analogously. 
Note that, in addition to (13.21) . F G and lim^-j-oo sup 3 ,g[g X)p|>d KP<z(x, f) = 0. 
Then, we deduce from m Theorem 8.2] (or [51 Theorem 11.2.3]) that 


N—^oo 


(3.8) 
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where V^° is the solution of the stochastic differential equation 

t 

r dF 

= Wc,it) - + J ^(^(^ 0 , v)) dv, (3.9) 

0 

with W, W being independent Brownian motions, and Ce(t) := J* q{z{zo,v), i)dv. 
Furthermore, V^° solves the following stochastic differential equation (see [131 Chap. 
9, Ex. 2.7]): 

t t 

= j \/ 9 {z{zo,v))dBy + J ^{z{zo,v)) VJ° dv, t>0. (3.10) 

0 0 

The last terms of the right hand sides of (13.91) and (|3.10ll have the same form. We 
claim that 

{Wc,it) - = ( j V9i<^o,v))dB 

Vo 

In this case, we deduce that == V^° (from the uniqueness of the solutions of 
(13.911 and (|3.10ll 'l. Hence, the first statement follows from (13.811 . In addition, since 
we have dlSl), and and z(zo,-) do not leave the interval [0,1], the remaining 
result is obtained from [in Theorem 8.5] (or [TUI Theorem 2.7]). 

Now we show the claim. The involved processes are local martingales with the 
same quadratic variation given by the deterministic function J* g{z{zo,v))dv. In 
addition, g{z{zo, v))dv = oo. The claim follows as an application of the Dambis- 
Dubins-Schwartz theorem (see for ex. [T31 Chap. 5, Theorem 1.6]). □ 



Lemma 3.5. The characteristic function ofV^° is given by 


/ 2 ( 20 ,^) \ 

exp |^-^F 2 (z( 2 :o,t)) / if zq ^ , 


= Xn ■ 


In particular, we have 

z{zo,t) 

Var(y/o) = <( I if^o^x+, 

.(i) |'l_,-2”vW*) ifzo = xl 


2VA 


(3.11) 


Proof. We give the proof for zq Xq , the remaining case follows similarly. The sec¬ 
ond statement follows from the first one, evaluating at 61 = 0 the second derivative 
of Ip with respect to 9. 

Note that the stochastic integral appearing in the definition of V^° is a local mar¬ 
tingale with deterministic quadratic variation given by f* g(z(zo, v))/F^(z(zo, v))dv. 
Since giz{zo, v))/F'^{z{zo, v))dv = 00 , we deduce from the Dambis-Dubins- 
Schwartz theorem (see [T31 Chap. 5, Theorem 1.6]) that 



y^g{z{zo,v)) 

F{z{zo,v)) 



t>0 



where (/3t)t>o is a standard Brownian motion. 


/ 


g(z(zQ,v)) 


dv 


0 ^ t>0 

The result is obtained using that 


□ 














THE DETERMINISTIC LIMIT OF THE MORAN MODEL 


7 


The main significance of Theorem (13.41) is that it provides a central limit theorem 
for the stationary distributions tt^w. We point out that such a result makes also 
part of nni Theorem 2.7] in the framework of density dependent families of Markov 
chains. However, in the latter there is a mistake in the variance of the corresponding 
limiting Gaussian distribution. For this reason, we provide a detailed proof to the 
next result. 


Corollary 3.6. Let be a random variable distributed as then 

In particular, we have ttzn —-—> 5^+. 

A-s-oo 


Proof. Note first that. Lemma 13.51 implies that limt_,.oo Var(V^^“) := In 

addition, the random variables are centred and Gaussian. We conclude that 




id) 

t—^oo 


' OO ’ 


(3.12) 


where Vf° 




addition, from the irreducibility of , we also have 




N 


id) 




N 


(3.13) 


where Z^ ~ tt^n. 

We fix now e > 0. Theorem 13.41 implies that, there is Nq^e) such that, for all 
N > iVo(£), supoQ \'ipN{t, d) — ipit, 0)1 < e. Therefore, we have for N > No(e) 


E 


- E 


< 

E 

giSVNZ^' 

- E 












^iS\/Nxf _ ^i0pNz{zo,t) 


E 


..levr 


- E 




Using (13.121) . (|3.13D . (I3.4|l and taking the limit when t tends to infinity in the 
previous inequality, we deduce that 


E 


„*eVlv(z^-x+) 


- E 




< e. 


Since this holds for all £ > 0, the result follows. 


□ 
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